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The original ideas about noncommuting coordinates are recalled. The connection between U(l) gauge fields 
defined on noncommuting coordinates and fluid mechanics is explained. Non-Abelian fluid mechanics is described. 



I The idea that configuration-space coordinates 
may not commute, 
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where 6 lJ is a constant, antisymmetric two- index 
object, has arisen recently from string theory, but 
in fact it has an older history. Like many inter- 
esting quantum-mechanical ideas, it was first sug- 
gested by Heisenberg, in the late 1930s, who rea- 
soned that coordinate noncommutativity would 
entail a coordinate uncertainty and could ame- 
liorate short-distance singularities, which beset 
quantum fields. He told his idea to Peierls, who 
eventually made use of it when analyzing elec- 
tronic systems in an external magnetic field, so 
strong that projection to the lowest Landau level 
is justified. But this phenomenological realiza- 
tion of Heisenberg's idea, which I shall discuss 
presently, did not address issues in fundamental 
science, so Peierls told Pauli about it, who in turn 
told Oppenheimer, who asked his student Snyder 
to work it out and this led to the first published 
paper on the subject pQ. 

The coordinate noncommutativity in the low- 
est Landau level is very similar to today's string- 
theory origins of noncommutativity - both rely on 
the presence of a strong background field. Also, 
thus far, it is the only physically realized example 
of noncommuting coordinates, so let me describe 
it in a little detail We consider the motion 

"This work is supported in part by funds provided by the 
US Department of Energy (DOE) under cooperative re- 
search agreement #DF-FC02-94ER40818 



of a charged (e) and massive (to) particle in a 
constant magnetic field (B) pointing along the 
z direction. All interesting physics is in the x-y 
plane. The Lagrangian for this planar motion is 

L = \m{x 2 + y 2 ) + e[xA x + yA y ) - V(x, y) (2) 

where the vector potential A can be chosen as 
(0,xB) and V(x,y) describes additional interac- 
tions ("impurities"). In the absence of V, the 
quantum spectrum consists of the well-known 
Landau levels |n, d), where n indexes the level's 
energy eigenvalue, and d describes the infinite de- 
generacy of each level. The separation between 
levels is 0(B/m), so that in the strong magnetic 
field limit only the lowest Landau level |0, d) is 
relevant. But observe that the large B limit cor- 
responds to small to, so projection on the lowest 
Landau level is also achieved by setting to to zero 
in J5J. In that limit the Lagrangian (J2J, in the 
chosen gauge, becomes 



LtLt = eBxy - V(x, y) 



(3) 



This is of the form pq — H{p,q), and immedi- 
ately identifies eBx and y as canonical conju- 
gates, leading in the usual way to the commutator 



\x, y] = —i—r 
1 ' yj eB 



(4) 



(The "Peierls substitution" consists of determin- 
ing the effect of the impurity by computing the 
eigenvalues of V(x,y), where x and y are non- 
commuting. There is much more to be said about 
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this. For example, the issue of ordering the (now) 
noncommuting arguments (x, y) of V must be set- 
tled. Also one needs to understand the behavior 
of wave functions. Before the (phase-space reduc- 
tive) strong- B limit they depend on both x and y, 
while after the limit they can depend only on one 
of the two noncommuting variables, all the time 
retaining their normalization. All this, as well as 
other matters, is explained in Ref. 

For another perspective, consider calculating 
the lowest Landau level matrix elements of the 
[x, y] commutator 

(0, d\xy - yx\0, d!) = M{d, d') - M*{d', d) (5) 
where 

M(d,d') = (0,d\xy\0,d') . (6) 

We evaluate (JBJ by inserting intermediate states 
in product xy: 

M(d,d') =^2(0,d\x\s)(s\y\0,d') . (7) 

s 

If the sum is over all the degenerate Landau lev- 
els, then one finds that vanishes: x and y do 
commute! But if one pretends that the world is 
restricted to the lowest Landau level and includes 
only that level (with its degeneracy) in the inter- 
mediate state sum 

M lLt = £)(Q, d|ar|0, d")(0, d"\y\0, d') (8) 

d" 

one finds that on this truncated state space, 
eq. © becomes consistent with |@J: 

(0,d\[x,y]\0,d') = -i^(0,d\0,d') . (9) 

A generalization now suggests itself: Suppose 
one retains not just the lowest Landau level, but 
includes the first N levels (setting N — for the 
lowest level) [3]. In this case one should consider 
all levels n < N, both as external states and con- 
tributing to the intermediate state sum in J7J. 
The result is that the matrix elements of the co- 
ordinate commutator vanish unless the external 
levels (n, n') coincide, n — n', and for the coinci- 
dent case the only nonvanishing matrix element 



is at the highest Landau level n = n' = N: 
(N,d\[x,y]\N,d') (10) 
= -i^(N + l)(N,d\N,d') . 

As N — > oo, more and more levels are included, 
the highest level decouples and the coordinates 
commute. 

II Let me now return to the general and abstract 
problem of noncommuting coordinates. When 
confronting the noncommutativity postulate l[T]l. 
it is natural to ask which (infinitesimal) coordi- 
nate transformations 

8x l = -f{x) (11) 

leave unchanged. The answer is that the (in- 
finitesimal) transformation vector function f l (x) 
must be determined by a scalar f(x) through the 
expression 0] 

f(x) = P*djf{x) . (12) 

Since then dif z (x) = 0, these are recognized as 
volume-preserving transformations. [They do not 
exhaust all volume-preserving transformations, 
except in two dimensions. In dimensions greater 
two, l|12|) defines a subgroup of volume-preserving 
transforms that also leave 8^ invariant.] 

The volume-preserving transformations form 
the link between noncommuting coordinates and 
fluid mechanics. Since the theory of fluid me- 
chanics is not widely known outside the circle of 
fluid mechanicians, let me put down some rele- 
vant facts 0. There are two, physically equiva- 
lent descriptions of fluid motion: One is the La- 
grange formulation, wherein the fluid elements 
are labeled, first by a discrete index n: X n (t) 
is the position as a function of time of the nth 
fluid clement. Then one passes to a continuous 
labeling variable n — > x : X n (t) — > X(t,x), 
and x may be taken to be the position of the 
fluid element at initial time X(Q,x) = x. This 
is a comoving description. Because labels can 
be arbitrarily rearranged, without affecting phys- 
ical content, the continuum description is invari- 
ant against volume-preserving transformations of 
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x, and in particular, it is invariant against the 
specific volume-preserving transformations 1|12|) . 
provided the fluid coordinate X transforms as a 
scalar: 

S f X = f\x)^-X = W'i),Xn,j . (13) 

The common invariance of Lagrange fluids and 
of noncommuting coordinates is a strong hint of 
a connection between the two. 

Formula 1)13(1 takes a very suggestive form when 
we rewrite it in terms of a bracket defined for 
functions of x by 

{dix), 2 {x)} = O^diO^dMx) . (14) 
Note that with this bracket we have 



(15) 



So we can think of bracket relations as classical 
precursors of commutators for a noncommutative 
field theory - the latter obtained from the former 
by replacing brackets by —i times commutators, 
a la Dirac. More specifically, we shall see that 
the noncommuting field theory that emerges from 
the Lagrange fluid is a noncommuting U(l) gauge 
theory. 

This happens when the following steps are 
taken. We define the evolving portion of X by 



X\t,x) = x l + e ij A j (t,x) 



(16) 



(It is assumed that 9^ has an inverse.) Then ((13J1 
is equivalent to the suggestive expression 



S f A i = a i f + {A i ,f} . 



(17) 



When the bracket is replaced by (— i) times the 
commutator, this is precisely the gauge transfor- 
mation for a noncommuting U(l) gauge potential 
Ai. Moreover, the gauge field fy emerges from 
the bracket of two Lagrange coordinates 

{X\ X J ) = 9 ij + 9 m 9 jn F mn (18) 
F m n = d m A n — d n A m + {A m , A n } . (19) 

Again (|19fl is recognized from the analogous for- 
mula in noncommuting gauge theory. 



Ill What can one learn from the parallel for- 
malism for a Lagrange fluid and a noncommut- 
ing gauge field? One result that has been ob- 
tained addresses the question of what is a gauge 
field's covariant response to a coordinate trans- 
formation. This question can be put already for 
commuting, non-Abelian gauge fields, where con- 
ventionally the response is given in terms of a Lie 
derivative Lf. 



8fx» = -F(x) 

S f A^ = L f A n = f a d a A u + d u f a A c 
But this implies 



(20) 
(21) 



SfF^ = LfF^ (22) 
= f ^ctF^y -\- Ofif F au -\- d v f F^ a 

which is not covariant since the derivative in the 
first term on the right is not the covariant one. 
Thus we consider the conventional approach to 
be defective because it does not preserve gauge 
invariance. The cure in this, commuting, situa- 
tion has been given some time ago |BJ: Observe 
that H21fl may be equivalently presented as 



SfA„ = L f A n = f a (d a A n - 8, L A a - i[A a ,A n }) 
+ f a d^A a - i[A», f a A a ] + d^f a A a 
= f a F a ^ + D^f a A a ) . (23) 

Thus, if the coordinate transformation generated 
by f a is supplemented by a gauge transforma- 
tion generated by —f a A a , the result is a gauge 
covariant coordinate transformation 



(24) 



and the modified response of F^ involves the 
gauge-covariant Lie derivative L'j\ 



5' f F nv = L' f F nv (25) 
— f F) a F nu -\- d^f F aif -\- dp f Ff ia . 

In the noncommuting situation, loss of covari- 
ance in the ordinary Lie derivative is even greater, 
because in general the coordinate transforma- 
tion functions f a do not commute with the fields 
A u , F nv ] moreover, multiplication of x-dependcnt 
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quantities is not a covariant operation. All these 
issues can be addressed and resolved by consid- 
ering them in the fluid mechanical context, at 
least, for linear and volume-preserving diffeomor- 
phisms. The analysis is technical and I refer you 
to the published papers |4I7| . The final result for 
the covariant coordinate transformation on the 
noncommuting gauge potential A^, generated by 
f a (x),is 

S' f \ = l{r(X),F afl } + (26) 
plus reordering terms. 

Note that the generating function f a (X) enters 
the anticommutator { , }+ with covariant argu- 
ment X. f a is restricted to be either linear or 
volume-preserving; in the latter case there are 
reordering terms, whose form is explicitly deter- 
mined by the fluid mechanical antecedent |4I7| . 

IV Next, I shall discuss the Seiberg-Witten 
map [Hj, which can be made very transparent 
by the fluid analogy. The Seiberg-Witten map 
replaces the noncommuting vector potential 
by a nonlocal function of a commuting poten- 
tial A^ and of 0] i.e., the former is viewed as a 
function of the latter. The relationship between 
the two follows from the requirement of stability 
against gauge transformations: a noncommuting 
gauge transformation of the noncommuting gauge 
potential should be equivalent to a commuting 
gauge transformation on the commuting vector 
potential on which the noncommuting potential 
depends. Formally: 

MA + dX) = A^ X \A) . (27) 

Here A is the Abelian gauge transformation func- 
tion that transforms the Abelian, commuting 
gauge potential A^; G(A, A) is the noncommuting 
gauge function that transforms the noncommut- 
ing gauge potential A^. G depends on A^ and A, 
and one can show that it is a noncommuting 1- 
co cycle [S]. 

Moreover, when the action and the equations 
of motion of the noncommuting theory are trans- 
formed into commuting variables, the dynamical 
content is preserved: the physics described by 



noncommuting variables is equivalently described 
by the commuting variables, albeit in a compli- 
cated, nonlocal fashion. 

The Seiberg-Witten map is intrinsically inter- 
esting in the unexpected equivalence that it es- 
tablishes. Moreover, it is practically useful for 
the following reason. It is difficult to extract 
gauge-invariant content from a noncommuting 
gauge theory because quantities constructed lo- 
cally from F^v are not gauge invariant. To 
achieve gauge invariance, one must integrate over 
space-time. Yet for physical analysis one wants 
local quantities: profiles of propagating waves, 
etc. Such local quantities can be extracted 
in a gauge-invariant manner from the physi- 
cally equivalent, Seiberg-Witten-mapped com- 
mutative gauge theory. 

For example, colleagues and I addressed the fol- 
lowing question Consider free, commuting 
Maxwell theory, and its noncommuting counter- 
part. The physical phenomena implied by the for- 
mer are well known: they are the electromagnetic 
waves, e.g., the plane waves. Can one find anal- 
ogous phenomena in the noncommuting theory, 
compare the two, and thereby set some experi- 
mental limits on the noncommutativity? 

The previously mentioned difficulty is now evi- 
dent: Waves are local excitations, but no meaning 
can be attached to local quantities in the non- 
commuting theory because they are gauge vari- 
ant. On the other hand, gauge-invariant noncom- 
muting quantities involve integration over space- 
time, but we would lose detailed physical infor- 
mation about Maxwell waves if their profiles are 
integrated. 

Here the Seiberg-Witten map suggests a solu- 
tion. Working to O(0), we express the noncom- 
muting fields in terms of commuting ones: 

A^ = A^-^A a (d p A tl -F^) (28a) 
% v = + B aSi F ail F 0v - 9 afs A a d p F^ (28b) 

where F^ is the Abelian field strength, the curl 
of A^. These profiles satisfy the condition (|27(1 
with 

G(A, A) = A - 9 afj A a df3\ . (29) 
Also <|28all solves to O(0) a differential equation 
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that can be derived from 1)27(1: 



i + 



OA - - - 



(30) 



Note however that when 9 a @ is larger than (2x2), 
(|3"U|) entails a set of equations, but they are not 
integrable. This means that when one constructs 
Antf ^ 0) from A^{9 = 0) = A^, the construc- 
tion depends on the path in 6*-space. We shall 
return to this point. 

The commutative action, which emerges when 
(|28bJ) is inserted in the noncommuting Maxwell 
action, reads 



(31) 



\9^F a pF^F^ - \e ap F m F u pF^ 



The "Maxwell" equations are 
d 



dt 



B + V x E = 
V -B = 



(32a) 
(32b) 



These are not dynamical; they merely express the 
fact that electric and magnetic fields {E % = F z0 , 
B l = —^e^ Fjk) are given in terms of potentials. 
The remaining equations, which follow from l(31() . 
can be given a Maxwell form 







— D - V x H = 

dt 



(33a) 

V • D = . (33b) 

with constituent relations following from l|31|) : 

D = (1 - 9 ■ B)E + {0 ■ E)B + (E ■ B)6 (34a) 
H = (1-0B)B 

+ \{E 2 - B 2 )6 - (9 ■ E)E . (34b) 

Here we have assumed that 9 a ^ has no temporal 
component: 9 0a = 0, = e ljk 9 k . 

It turns out that [to 0{9)] the nonlinear equa- 
tions (|32"|l - il34*|) still possess plane wave solu- 
tions. However, in the presence of a background 



magnetic field b the dispersion law between fre- 
quency lu and wave number k is modified from 
the Maxwell result u — k; we find 



u) = k(l - 9t ■ b T ) 



(35) 



where the subscript "T" denotes components 
transverse to the direction of propagation k. 

Evidently, Lorentz symmetry is violated: the 
velocity of light has changed, increasing or de- 
creasing over the Maxwell value. In principle 
this should lead to a fringe shift in a Michaelson- 
Morley-type experiment. However, with present 
limits on 9 < (10 TeV)~ 2 (obtained from stud- 
ies of interaction with matter ^2] ), visible light 
(wavelength ~ 10~ 5 cm) and attainable magnetic 
fields (~ 1 tesla) would lead to the shift of 1 fringe 
on an optical length of 1 parsec (10 18 cm). Ob- 
viously more advanced technology is required to 
render feasible this test of noncommutativity! 

Note that both polarizations undergo the same 
velocity change. This is in contrast to a previous 
model for violating Lorentz symmetry in electro- 
dynamics, wherein a spatial Chern-Simons was 
added to the Maxwell action [2]. In this modifi- 
cation the theory remains linear; eqs. ((33(1 involve 
E and B [i.e., the constituent relations l(34|l take 
on their 9 = value] but the right side of (|33a|l 
acquires a term proportional to B. The net re- 
sult is that plane waves with different polariza- 
tions travel with different velocities, producing a 
Faraday-type rotation, which can be ruled out ex- 
perimentally |13I14| . 

Although a Faraday rotation does not occur 
in the noncommutative generalization of Maxwell 
theory, 1)31 [I. the following is noteworthy. Ob- 
serve that the modification in 1)31(1 involves two 
numbers (1/8 and 1/2). Since 9 may be rescaled, 
only their ratio (@J matters. One can check that 
when the ratio departs from this particular value, 
plane waves still solve the equations, but now the 
two polarizations travel at different velocities. So 
noncommutativity is the unique modification to 
Maxwell theory [within the class l(31(l ] that avoids 
a Faraday rotation. 

A final comment: It has been observed that 
the action 1(31(1 may be viewed as a Maxwell ac- 
tion in a background gravitational field, which it- 
self is constructed from F^ v . Moreover, the equa- 
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tions l|32[l - i|34[) coincide with the geodetic equa- 
tion in that gravitational field |15| . This hints at a 
view that noncommutativity reflects the presence 
of a "medium" , and this theme will reappear in 
an explicit realization of the Seiberg-Witten map 
within fluid mechanics. 

V The next investigation I shall describe re- 
sults in an explicit formula for the Seiberg-Witten 
map. The derivation again makes use of the fluid 
analogy, but now we need the second, alterna- 
tive formulation of fluid mechanics, the so-called 
Euler formulation. This is not a comoving de- 
scription (like the Lagrange formulation), rather 
the experimenter observes the fluid density p and 
velocity v at a given point in space-time (t,r). 

In preparation for our derivation, I shall briefly 
outline the Euler point of view. Moreover, later I 
shall rely on that outline to give a generalization 
of Eulerian fluid mechanics to the case when a 
non-Abelian symmetry group acts on the various 
degrees of freedom (components) of the fluid. 

A point particle in space-time is described by 
its space-time coordinat X M (r), which depends 
on an evolution variable t that parameterizes the 
path. One expects that X^ satisfies some dynam- 
ical equation, which determines X^, but rather 
than concentrating on that equation (this would 
correspond to a Lagrange description) we focus 
on the density and current associated with that 
particle. These are functions of x^ = (t, r) and 
are related to X^ by the formulas 



p(t,r) = J dTX°(T)S 4 (X^{T) -x") (36a) 
= J dTX°{T)S(X (T)-t) S 3 (X(T)-r) 

j(t,r)= J dTX(T)S 4 (X^(T) -a") (36b) 

= J dTX(r)S(X°(T)-t) S 3 (X(T)-r) 

We can choose the parameterization so that 
X°(t) — r, whereupon l|36l) becomes 



p(t,r)=6 3 (X(t)-r) 
j(t,r)=X(t)6 3 (X(t)-r) 



(37a) 
(37b) 



Observe that either from the definition H36JI or 
from (|3T| a continuity equation for (p,j) follows 
identically: 

^- t p(t,r) + V -j(t,r) = . (38) 

Note also that l)37bjl defines the velocity v as a 
function of (t, r): it is X(t) with X(t) replaced 
by r: j = pv. 

A second equation is obtained by differentiat- 
ing i|37b|l with respect to time: 

d ■ • ■ • __ , 

j = pv + pv = pv — vv ■ (pv) 

= XX 3 -^-S 3 (X{t) -r)+ force . (39a) 



Of 



Here "force" denotes a formula for X , which spec- 
ifies the dynamics that govern particle motion. 
The next to last term may be rewritten as 

-XX J -^-8 3 (X(t) -r)= -^(vv j p) . 
Thus putting everything together we find 



d 1 r 

— v + v ■ vd = - force. 
at p 



(39b) 



Eqs. H38fl and (|39 b|) are the Eulerian equations 
for a "fluid" composed of a single particle! When 
there are several particles, labeled by index n, 
formulas (|37fl are replaced by 



p(t,r)=Y, 63 ( X ^- r ) 

n 

j{t,r) = p(t,r)v(t,r) 

= S 3 (X n (t)-r) 



(40a) 



(40b) 



but the subsequent development is as before, cul- 
minating in and (|39b(l . Finally, for a fluid, 
the discrete label becomes a continuous parame- 
ter specifying the fluid n — > x, X n (t) — > X(t,x); 
the density and velocity expressions are the con- 
tinuum generalizations of l|40|) 

p(t,r) = J d 3 x5 3 (X(t 7 x) -r) (41a) 
j(t,r) = p{t,r)v(t,r) 

d 3 xX(t,x)8 3 (X(t,x) -r) (41b) 



7 



and once again the continuity l|38fl an d Euler 
(|39b(l equations are established. 

It is useful to understand in detail the content 
of (|41|l , which expresses the relation between the 
Lagrange variables and the Euler variables 

[p(t,r),v(t,r)]. Evidently the integration over 
the parameter x evaluates the delta function at a 
value of x = \(t, r) such that 



(42) 



and also there is a Jacobian. Thus (|41a|l also 
states that 



1 



det 



dX%x) 



p(t, r) dxi 
while (|41b(l implies 



x=x 



v(t,r) = X\t,x) 



(43) 



(44) 



Effectively the passage from Lagrange to Euler 
descriptions involves interchanging the dependent 
variable X with the independent parameter x, 
which is renamed r with the help of x{t, r )- 

Finally we note that the theory is completed 
by giving a model for the force. This will include 
the (negative) gradient of the pressure, frequently 
taken as a function of p. If the fluid is charged, 
we are dealing with megnetohydrodynamics: 1/p 
force includes the Lorentz force: E + v x B, and 
j 1 * = (p, pv) is the source for the electromagnetic 
fields in Maxwell's equations 



(45) 



Later I shall describe how this is generalized when 
an internal symmetry is present, and a group in- 
dex a labels various physical quantities. 

Both a configuration space Lagrangian/action 
formulation and a canonical Hamiltonian formu- 
lation for equations 1(38(1 , (|39bjl with pressure and 
Lorentz forces, and 1)45(1 can be given 0. 

VI Let me now turn to the task of obtaining 
an explicit formula for the Seiberg-Witten map 
from the fluid analogy. Actually, I shall present 
the inverse map, expressing commuting fields in 
terms of noncommuting ones; also the derivation 



will be restricted to the case of two spatial di- 
mensions [where 6 lJ involves a single quantity and 
there are no integrability conditions on the differ- 
ential equation l(3U(l ]. Higher dimensional cases 
have been been treated as well; the development 
is more complicated. A few remarks about them 
will be given later; for a complete discussion I 
refer you to the published literature |4j. 

The (inverse) Seiberg-Witten map, for the case 
of two spatial dimensions, can be extracted from 
H41[) . which we rewrite as 

j»(t,r) = J d 2 x^- t X»(t,x)5 2 (X(t,x)-r) 

[X° = 1] . (46) 

Observe that the right side of (|46|l depends on 
A through X [sec l|16f) ]. It is easy to check that 
the integral 1(461) is invariant under the transfor- 
mations (|13|) : equivalently viewed as a function 
of A, it is gauge invariant [see ifTTj) ]. Owing to 
the conservation of j^c^j 7 ' = [see its 
dual Eq/^j 71 satisfies a conventional, commuting 
Bianchi identity, and therefore can be written as 
the curl of an Abelian vector potential A a , apart 
from proportionality and additive constants: 



d a A/3 — dpAa + constant 

oce^y d 2 x^- t X»5 2 (X-r) 
diAj — djAi + constant 

J d 2 xS 2 (X -r) = E i:j p 



(47) 



This is the (inverse) Seiberg-Witten map, relating 
the A to A. 

Thus far the operator properties and the non- 
commutativity have not been taken into account. 
To do so, first we reinterpret the integral over 
a; as a trace over the Hilbert space that realizes 
the noncommutativity of the coordinates (this is 
a standard procedure); second we must provide 
an ordering for the (5-function depending on the 
operator X' 1 = x l + O l ^Aj. This we do with the 
Weyl prescription by Fourier transforming. The 
final operator version of equation l|47|) . restricted 
to the two-dimensional spatial components, reads 
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A 2 re lk r {d iaj - djca) 



d 2 xe lk x - (2ir) z 5{k) 



(48) 



Here the additive and proportionality constants 
are determined by requiring agreement for weak 
noncommuting fields. 

Formula l|48l) has previously appeared in a di- 
rect analysis of the Seiberg-Witten relation |16|. 
Now we recognize it as the (quantized) expression 
relating Lagrange and Euler formulations for fluid 
mechanics 

In higher-dimensional cases, greater than 2, 
there arises the need to define a generalized cur- 
rent tensor, whose dual gives rise to a two-index 
tensor satisfying a Bianchi identity (closed and 
exact 2-form). There is some ambiguity in the 
construction when quantum mechanical ordering 
needs to be implemented. The ambiguity pre- 
cisely parallels the ambiguity in "integrating" the 
path-dependent Seiberg-Witten differential equa- 
tion; see Ref. @|. 

VII A natural question that comes at this 
stage is whether one can extend the fluids- 
noncommuting fields analogy to include a non- 
Abelian symmetry group. This has not been 
achieved thus far, but for a first step we need 
to construct non-Abelian fluid mechanics. This 
has been done. I shall now present the the- 
ory, which arises very naturally from non-Abelian 
point-particle mechanics, just in the same was as 
I presented earlier a derivation of an (Eulerian) 
fluid from point-particle mechanics |17| . 

Clearly, a non-Abelian fluid will possess a cur- 
rent J£ that carries an internal symmetry index a. 
In the presence of dynamical gauge fields (non- 
Abelian magnetohydrodynamics) J£ serves as 
source for these fields in a generalization of (|45|l . 



{D^)a = J V a 



(49) 



and is covariantly conserved, for consistency with 



Here is the covariant derivative (D^) ac = 
d fl S ac + f A° , where the non-Abelian gauge po- 
tential Af, is contracted with the structure con- 
stants f abc of the non-Abelian group. The chal- 
lenge for us is to discover how the non-Abelian 
current determines a charge density and velocity, 
and what kind of Euler equations are satisfied by 
these quantities. 

To discover this, we recall the dynamics of 
a non-Abelian point particle, as described by 
Wong PH]- The particle carries a coordinate X(t) 
and a non-Abelian charge q a (t)- The charge and 
current densities read [compare i|37|) ] 



p a (t,r)=q a (t)6 3 (X(t)-r) 
j a (t 7 r)=q a (t)X(t)5 3 (X(t) 



(51) 



Covariant conservation is assured because the 
charge q a {t) is assumed to obey the Wong equa- 
tion 



q a (t) + r bc X»(t)A b ^(t,X(t))q c (t) = 
[A = 1] . 



(52) 



Moreover, particle acceleration is determined by 
the equation satisfied by X(t), e.g., 



Pfj, — qaF^x" 



(53) 



( V)a = 



(50) 



where is the particle 4-momentum. A La- 
grangian/action and a Hamiltonian canonical for- 
mulation of these equations can be given; it in- 
volves the Kirilov-Kostant 1-form on the Lie al- 
gebra [T§] . 

The passage to a fluid is now clear. First 
we consider a finite number of particles, and 
introduce a particle label n: X^,P",g™. 
Next we pass to the continuum n — * x: 
X^(t, x), Pfi{t, x), q a (t, x) - these are the La- 
grange variables. The non-Abelian charge and 
current densities are defined as in (|41|l 



Pa (t,r)= I d 3 xq a (t,x)6 3 (X(t,x)-r) (54a) 
3 a (t,r)= (54b) 
d 3 xq a (t,x)X(t,x)S 3 (X(t,x) -r) . 
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This remains covariantly conserved provided the 
Lagrangian non-Abelian charge variable q a (t,x) 
satisfies 



d_ 

Of 



q a (t,x) 



(55) 



+ f abc X»(t, x)A b fl (t, X(t, x))q c (t, x) = 



The important feature of these formulas is that 
by virtue of the <5-function they factorize into an 
Abelian current and a non-Abelian charge factor 



Pa(t,r) = Q a (t,r)p(t,r) 
3 a (t,r) = Q a (t,r)j(t,r) 

= Qa(t,r)p(t,r)v(t,r) 



(56a) 



(56b) 



with Q a being the Eulerian charge, related to the 
Lagrange charge by 



Q a (t,r) = q a (t,x) 



=x(t-r) 



(57) 



[see (|42I) - H44[) ]. Moreover, the Abelian current is 
still given by (|41|1 . so it satisfies the Abelian con- 
tinuity equation, and v satisfies an Euler equa- 
tion that reflects the underlying particle dynam- 
ics. Finally, the Eulerian charge (|56l) obeys the 
fluid Wong equation that follows from l|55|) and 
()57J) or alternatively from (|56J) when it is remem- 
bered that j£ = Qaj^ and j M are respectively 
conserved covariantly and ordinarily: 



d_ 
di 



Qa+V- VQ a = fabcQ b (A C - V ■ A C ) 



(58) 



Non-Abelian magnetohydrodynamics then is de- 
scribed by the field equations 149fl . I|50(l . by the 
factorization l|56|l . by the fluid Wong equation 
{SB} ; an d by the Euler equation 

d 1 
— v + v ■ Vv = q a E a + v x q a B a + -force (59) 
ot p 

where we have explicitly exhibited the contribu- 
tion to the force term of the non-Abelian Lorentz 
force, involving the non-Abelian electric E a and 
magnetic B a fields. A Lagrangian for these equa- 
tions makes use of a field-theoretic generalization 
of the Kirilov-Kostant 1-form |17|. 



In the above approach a single density and ve- 
locity describe all the group components, i.e., p 
and v are "a" -independent. It is possible to give 
a more elaborate treatment where different p and 
v are associated with different elements of the 
Cartan subalgebra of the group. Thus for SU(2), 
with one Cartan element, there would still be only 
a single density and velocity describing all three 
group degrees of freedom. On the other hand, for 
SU(3) the more refined treatment uses two dis- 
tinct densities and velocities. For this elaboration 
I refer you to the literature |17| . 

Our approach to fluid was motivated by an un- 
derlying particle picture, as in the passage from 
<|37|l to (|4*T1) . There is also a field-based approach, 
which is illustrated by the fluid interpretation of 
the Schrodinger equation, due to Madelung |2()j . 
Observe that the Schrodinger equation for a par- 
ticle in an electromagnetic field derived from po- 
tentials (<p, A) reads 



eA) 2 ^ 



(60) 



If we define ^ = p l / 2 e l6 / h , then the real and 
imaginary parts of l|6UI) become 



d_ 

di 

d 

_,Vp ! (Vp) 2 



p + V-((V9-eA)p) 

Uve-eA) 2 



-ecp + % M 



M 



P 



8 p 2 



(61a) 
(61b) 
(61c) 



The first of these suggests identifying V6> — eA as 
the velocity, so that the continuity equation 
is regained. Then upon taking the gradient of the 
second equation, we arrive at 

d 1 
— v + v • = e(E + v x B) + -force (62) 

where "force" is the so-called Madelung 0(h 2 ) 
quantum force given by pVM. Thus we see that 
Eulerian fluid mechanics emerges within this re- 
analysis of Schrodinger theory. 

One may therefore assay a field-based non- 
Abelian fluid mechanics by beginning with a 
Schrodinger theory for a multicomponent wave 
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function transforming under some non-Abelian 
group. What emerges is a fluid mechanics that is 
much more complicated than the particle-based 
model that we presented above. The crucial dif- 
ference is that here there is no simple factoriza- 
tion of the current, nor does a Wong equation 

hold HMD- 

It is unclear at present which (if any) of these 
non-Abelian generalizations for fluid mechanics 
will prove to be the most useful. 
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